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When a crystalline solid such as iron is cooled 
across a structural transition, its final microstruc- 
ture depends sensitively on the cooling rate H3. 
For instance, an adiabatic cooling across the tran- 
sition results in an equilibrium 'ferrite', while a 
rapid cooling gives rise to a metastable twinned 
'martensite' [g]. There exists no theoretical 
framework to understand the dynamics and con- 
ditions under which both these microstructures 
obtain. Existing theories of martensite dynam- 
ics ®^$ describe this transformation in terms of 
elastic strain, without any explanation for the oc- 
curence of the ferrite. Here we provide evidence 
for the crucial role played by non-elastic variables, 
viz., dynamically generated interfacial defects. A 
molecular dynamics (MD) simulation of a model 
2-dimensional (2d) solid-state transformation re- 
veals two distinct modes of nucleation depending 
on the temperature of quench. At high temper- 
atures, defects generated at the nucleation front 
relax quickly giving rise to an isotropically grow- 
ing 'ferrite'. At low temperatures, the defects 
relax extremely slowly, forcing a coordinated mo- 
tion of atoms along specific directions. This re- 
sults in a twinned critical nucleus which grows 
rapidly at speeds comparable to that of sound. 
Based on our MD results, we propose a solid- 
state nucleation theory involving the elastic strain 
and non-elastic defects |fr|g|], which successfully 
describes the transformation to both a ferrite and 
a martensite. Our work provides useful insights 
on how to formulate a general dynamics of solid 
state transformations. 

We study a 2d model system which exhibits exactly 
two distinct equilibrium solid phases — a square and a 
triangular phase (with a rhombic unit cell) using an MD 
simulation on particles ('atoms') interacting with an ef- 
fective 2- and 3-body potential (see Fig. 1). By tuning the 
potential parameters, we induce a square-to-triangular 
solid-state transformation at constant temperature. Fig- 
ure 1 shows two typical 'quenches' (referred to as high 
temperature T\ and low temperature T 2 quenches) start- 
ing from the equilibrium square phase. In both cases 
following the quench, the product phase is formed by 
a process of nucleation and growth (Fig. 1 (inset)), but 
while the T\ nucleation is homogeneous, the T 2 nucleation 
needs to be seeded by a single defect (heterogeneous). 
Nucleation of a triangular region within the square phase 



gives rise to vacancies/interstitials due to atomic mis- 
match, which relax with a temperature dependent dif- 
fusion coefficient W. Note that (a) the use of effective 
potentials rather than 'realistic' atomic potentials allows 
us to probe late times and large sizes ; (b) the experimen- 
tally more accurate quench schedule is a cut through the 
T—V3 plane, however the essential physics is unaltered by 
our constant T quench allowing for a more transparent 
interpretation of the MD results. 

In an earlier paper ||, we had presented a coarse- 
graining procedure which starts from a microscopic de- 
scription of a solid and constructs continuum variables 
relevant for a description of the solid at larger length 
scales. Thus we define coarse-grained dynamical elastic 
variables from the instantaneous displacements u of the 
particles from the ideal square lattice 11 - - nonlinear 
strain ey = l/2(diUj + djUi + diUidjUi), density fluctua- 
tion about the mean (p), c/> = (p— (p))/{p) (we shall refer 
to this as the vacancy/interstitial field) and dislocation 
density (density of 5-7 disclination pairs). 

In the T\ quench the vacancies/interstitials diffuse fast 
and annihilate each other, thus encouraging the trian- 
gular phase to grow isotropically (Fig. 2a) with a size 
R that grows linearly with time t. The critical nucleus 
is untwinned leading to a polycrystalline triangular solid 
('ferrite'). 

In the T 2 quench the seed vacancy does not diffuse 
away ; dominant fluctuations correspond to a coordinated 
movement of a line segment of atoms along the square 
axes (Fig. 3a). Once these fluctuations get large enough, 
they grow into a twinned critical nucleus with the twin- 
ning plane along one of the square axes. The twinned 
nucleus grows anisotropically, with a constant tip speed 
roughly equal to that of longitudinal sound (Fig. 3b). 
The vacancies and interstitials while staying separated 
from each other, are found concentrated at the interface 
of the growing nucleus (Figs. 3c). In the direction trans- 
verse to the twinning 'plane', the nucleus grows by step 
formation (Fig. 3a). 

The MD simulation highlights the crucial role played 
by defects in determining the mode of nucleation and the 
nature of the critical nucleus. From these simulations we 
learn that the relevant variables to describe the dynam- 
ics of nucleation should include both the strain ey pQ] 
and the defect density 4>. We now describe an analytic 
theory of solid-state nucleation for which the simplest 
dimensionless free-energy functional describing the first- 



order square (< e xy >= 0) to rhombus 
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transition at a fixed degree of undercooling a has the 
form, 
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where Wij = diUj is the deformation tensor. A quench 
across this transition to where the square phase is 
metastable, nucleates a small 'droplet' of the rhombus. 
Following the MD, we simultaneously create an envelope 
of 4> around the droplet ; this is the initial condition for 
the dynamical equations in the displacements u Jf|Jll[ 
and <f>: d t <t> + v • Vtf> = D v W 2 (2jcf> + Ccu), where the mi- 
croscopic vacancy diffusion coefficient has an Arrenhius 
form D v = Doo exp(— A/ksT), with an activation energy 
A §. The velocity field v = li + TCV0, where T is the 
microscopic kinetic coefficient of u. 

For a high temperature quench D v is large, so <j) re- 
laxes much faster than the time for the critical nucleus to 
form ; subsequently the nucleus grows isotropically into 
an equilibrium triangular solid ('ferrite') MM with R~t. 

For intermediate temperature quenches D v is smaller ; 
the nature of the critical nucleus is determined by com- 
puting the time taken for critical nucleation r„ ('first- 
passage-time'). Apriori we do not know whether the 
critical nucleus is twinned or not, and so we perform a 
variational calculation for both cases ; the 'true critical 
nucleus' is the one for which r n is smaller. 

In the linearised equation for </>, the velocity field (ig- 
noring lattice vibrations) reduces to the front velocity 
v = VfhS(r — R(£)), where R(i) is the position of the 
moving interface and n is the unit outward normal to 
the interface. Thus in the reference frame of the inter- 
face, the dynamics of <f> is diffusive. Since the relaxation 
time to the local minima of Eq. (pi) is smaller than bar- 
rier hopping times, we may treat e^- as 'slaved' to <fi 
in the calculation of the barrier height |12|. We then 
use a Kramers' formula to evaluate r„, thus apart from 
unknown prefactors, r„ = T _1 exp(A_E* /fcsT), where 
AE* is the </>(r, independent critical barrier energy. Fig- 
ure 4 clearly shows that for high temperature quenches 
(low undercooling), the true critical nucleus is a 'ferrite', 
while for a low temperature quench it forms a twinned 
bicrystal. We successfully reproduce two distinctive fea- 
tures of the well known 0% isothermal transformation 
curves of martensites — the horizontal transformation 
curve beyond a well defined 'start' temperature M s in- 
dependent of Doo, and a ferrite nose. It should be re- 
marked that such a calculation could not even be ad- 
dressed within the conventional viewpoint. 

What happens after the twinned critical nucleus is 
formed ? It is difficult to say based on our early time 
MD studies, whether the growing nucleus will eventu- 
ally add on more twins or whether it will continue as 
a bicrystal with a single twin interface. On the other 
hand we may analyse the late time dynamics within our 



analytic formalism. In the limit D v — ► when <fi re- 
mains frozen at the moving interface, a constrained vari- 
ational calculation |?]]|] , shows that the growing nucleus 
finds it favourable to add more twins. The free en- 
ergy of an anisotropic inclusion of length L and width 
W, with N twins along the length is given by F = 
AFLW + a pp (L + W) + a tw NW + [3(L/N) 2 N, where 
AF < is the free energy difference between the square 
and triangular phases, o~ pp and o~t w are the surface ten- 
sions of the square-rhombus and the twin interfaces re- 
spectively, and the last term is the <j> contribution at 
the square-rhombus interface. Minimization with respect 
to N gives L/N ~ W 1 ^ 2 , a relation that is empirically 
known for 2d martensites like In-Tl Binl . Moreover the 
energy of the nucleus is smaller when oriented along spe- 
cific directions, the habit 'planes' ||. 

The square-rhombus interface is studded with an ar- 
ray of vacancies such that </> —> on an average, with a 
separation equal to the twin size. Thus the strain decays 
exponentially from the interface over a distance of order 
L/N. This 'fringing field' [Q arises dynamically in our 
calculation, rather than from an imposed boundary con- 
dition pp[ . Indeed throughout our analysis we do not 
impose specific interfacial conditions PJlO|| , preferring to 
allow the solid to choose its own dynamical path. 

The theoretical formalism just presented may be gener- 
alised to aribitrary solid-state transformations. We have 
found that in order to develop a general theory for solid 
state nucleation, it is necessary to go beyond nonlinear 
elasticity theory. More specifically, the complete set of 
relevant dynamical variables should include non- elastic 
degrees of freedom like defect fields. Such non-elastic 
degrees of freedom should be relevant whenever a solid 
undergoes large shear deformations. Our work provides 
a framework to understand the dynamics of solid state 
transformations and the conditions under which various 
microstructures obtain. The methodology of translating 
from microscopic simulations to coarse-grained dynami- 
cal fields should be useful in other problems in material 
science such as brittle-ductile transitions. We are cur- 
rently generalising our theory to include the dynamics 
of dislocations, substituitional and interstitial impurities 
and a coupling to external stress to study the dynamics 
of shape-memory alloys. 
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Figure Caption 

Fig. 1 Phase diagram in the T — v$ plane at a 

density (p) = 1.05 using an MD simulation in the NVT 
ensemble. The particles interact via an effective poten- 
tial consisting of a 2-body V2 {rij ) = «2 ( fJ / r ij ) and 
a short range 3-body V^r^r.,-, r fe ) = v 3 [sin 2 (4% fc ) + 
sin 2 (49 jki) + sin 2 (40fcy)]. Particles i and j are separated 
by a distance ry. 9ijk is the bond angle at j between 
triplets (ijk). V2 favours a triangular lattice ground state 
while V3 is minimised when 9ijk = 0, 7r/4, tt/2 and so 
favours a square lattice. The units of length and energy 
are set by a and v-i respectively, making the unit of time 
a^/m/v2, where m is the particle mass. The MD time 
step is 0.001 corresponding roughly to a real time of 1 fs. 
The order parameter (tt) =(fi N)- 1 X^ fc ( sin2 ( 4 %fc)) 
takes values 0(1) in the square (triangular) phases respec- 
tively. The phase diagram was obtained by equilibrating 
a system of N = 1024 particles at various values of T and 
V3. Arrows marked T\ and Ti are constant temperature 
quenches at T\ = 1 and T 2 = 0.15. The values of u 3 to 
which the solid is quenched is W3 = 3 and v 3 — 1.2 respec- 
tively, which lie in the region of the phase diagram where 
the square is metastable (within the dotted line). Inset 
shows the time evolution of (f2) for the two quenches, 
T\ (thin line) and T2 (thick line), indicating a typical 
nucleation process. 

Fig. 2 MD snapshot at t = 10 4 following a T\ 
quench: 'Ferrite' nucleus grows isotropically as a poly- 
crystalline triangular region. The colours code the local 
values of (O) going from red in the square region to blue 
in the triangular region. Inset shows x-y trajectories for 
5 x 10 3 < t < 25 x 10 3 of 5 particles chosen to lie along 
a row (x-axis) at t = 0. At these times the particles are 
part of the growing nucleus, and their trajectories reveal 
significant diffusive motion. Number of unit cells of the 
square in each direction is 110. The vacancy/interstitial 
density <j) relaxes fast and quickly averages to 0. 

Fig. 3 MD snapshot at t = 6 X 10 3 following a 
T2 quench: (a) Twinned nucleus grows anisotropically. 
Colour coding as in Fig. 2. Inset shows trajectories for 
10 3 < t < 16.5 x 10 3 of 5 particles chosen to lie along 
a row (jc-axis) at t = (movement from bottom to top 
of inset). At these times the particles are part of the 
growing nucleus, and their trajectories reveal highly co- 
ordinated, 'military' motion characteristic of martensites. 
Number of cells is the same as above, (b) Positions of 
5 particles chosen to lie along a twin-column (j/-axis) as 
a function of time. Apart from thermal oscillations, the 
particle postions change only when they encounter the 
moving front at t = t s . A plot of the position of the 
ith particle yi as a function of t s gives the front-velocity 
Vf = 9.7. This compares well with the sound velocity 
calculated from the bulk modulus in the square phase, 
(c) The vacancy/interstitial density <p integrated over x 
as a function of y (along the twinning plane) at times 
t\ = 2000 and £2 = 5000. The arrow shows the direction 
of the displacement field resulting in an increased density 



at one end of the nucleus and a reduced density at the 
other. Inset shows the shear strain e xy integrated over y 
showing the twin boundary at ti, fa- 
Fig. 4 First-Passage-Times versus degree of un- 
dercooling a for a twinned and untwinned critical nu- 
clei. The calculation proceeds by recognising that a per- 
fect triangular solid is obtained from a perfect square by 
the deformation, R' t = (&y + eij)Rj, involving a shear 
f-xy = tyx = e and volume compression e xx — e yy = e 2 /2. 
We therefore parametrize e^ by the single function e(r, t). 
This admits a variational ansatz for (1) 'ferrite': e = 
outide a grain of size L and e = eo inside, with a smooth 
interpolation of width £ and (2) 'twinned nucleus': e = 
outide a rectangular grain of length L, width W and 
e = ±eo inside, connected by smooth interpolations. The 
diffusion equation, dt<p = D v V 2 (2-f(j> + Ceu) is solved 
with the initial condition cj>(r, 0) = Au ■ n (Au is geo- 
metrical mismatch at the square-triangle interface and is 
proportional to the length of the interface). Note that 
Au = at the twin interface. The free-energy of the 
grain at time t, E(L, W, t ; {&}) is obtained from Eq. (0), 
which we minimize with respect to the widths {£,} for 
every L, W and t. The barrier energy AE* and size 
L*,W* of the critical nucleus at every time t, is deter- 
mined by the saddle-point. The energy of the critical 
nucleus decreases with time, and so a crude but easily 
calculable estimate of the first-passage-time is obtained 
by a self consistent solution of r„ = T _1 cxp(AS*(r„)) . 
The curves have been calculated for Doo = 10 14 , A = 7, 
7 = 0.2 and C = 0. The upper dotted line represents the 
equilibrium transition temperature (zero undercooling). 
At small undercooling <\> relaxes fast and the critical nu- 
cleus is a 'ferrite', the first-passage-time for which goes 
through a minimum. At larger undercooling, the sin- 
gle twin nucleates faster (bold line) than the untwinned 
grain (lower dotted line) . The arrow shows the degree of 
undercooling M s below which the twinned nucleus forms 
('martensite start' temperature). 
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